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Abstract - The quality of the electrical contact between a 
cultured neuron and a substrate embedded 
microelectrode is of importance for effective transfer of 
an extracellular applied stimulus current to the 
intracellular potential. It is affected by the resistance of 
the seal, i.e. the gap between the cell membrane and the 
substrate, which restricts the leakage current, thereby 
favouring the efficiency of the stimulation current. 
The effects of variations in the geometry of the neuron- 
electrode interface on the sealing resistance and on the 
stimulus transfer are studied using a finite element model 
of this interface. Variations in the geometry of the 
neuron-electrode interface are represented by the 
excentricity, x c , of a pillbox shaped neuron with radius r„ 
cultured on an electrode with radius r e . 
The results indicate a sharp decrease in both sealing 
resistance and stimulus transfer when a transition occurs 
from complete sealing to defect sealing. At that point the 
leakage current splits up into a current through the gap 
and a current through the sealing defect. 
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I. Introduction 

Planar substrate electrode arrays (PEAs) offer the possibility 
to contact neuronal cultures with a large number of cell-size 
electrodes [1][2]. In our project, PEAs are used as a tool for 
the development of a cultured neuron probe, a new type of 
neuro-electronic interface for highly selective stimulation 
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purposes. 

For effective transfer of an extracellular applied stimulus 
current to the intracellular potential, the quality of the 
electrical contact between a cultured neuron and an electrode 
is of importance [3]. This quality is affected by the resistance 
of the seal, i.e. the gap between the cell membrane and the 
substrate, which restricts the leakage current, thereby 
favouring the efficiency of the stimulation current (fig. la). 
In previous work, a method for measurement of the sealing 
resistance in the neuron-electrode interface was proposed as a 
first step towards assessment of the quality of this contact [4], 
This method was based on a model of a neuron-electrode 
interface geometry in which the electrode was completely 
covered, i.e. sealed by the neuron. 

Now, this model is extended to include the effects of varying 
neuron-electrode geometries on the transfer of extracellular 
applied stimuli. Variations in the geometry of the neuron- 
electrode interface are represented as an excentricity, x c , of a 
neuron with radius r c , cultured on an electrode with radius r e 
(fig. lb). When x c increases from zero, the sealing resistance 
will decrease due to the altering geometry of the gap between 
the cell membrane and the substrate. At x c =r c -r e , a transition 
occurs from complete sealing to defect sealing: The leakage 
current splits up into a current through the gap and a current 
through the sealing defect. Further increase of x c will also 
decrease the area of the patch membrane (fig. lc). 
The effects of the neuron excentricity on the sealing resistance 
and on the stimulus transfer are computed using finite element 
analysis. The transitions are studied at several geometries by 
varying the radii of neuron and electrode. The results indicate 
a sharp transition in both sealing resistance and stimulus 
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fig. 1 (a) Schematic presentation of the neuron-electrode interface. The efficiency of the stimulation current is increased by restriction oj 
the leakage current through the. gap. (h) Model assumption on the geometry of the neuron-electrode interface: the neuron is positioned 
excentrically on the electrode (with excentricity x c ). (cj Electrical equivalent circuit of the (passive) neuron-electrode interface. The sealing 
resistance and area of the patch membrane are modulated by the excentricity, x„ of the neuron. 
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transfer around x,;=r c -r e . In case of complete sealing the 
stimulus transfer remains relatively constant. 



II. Methods 

In fig. 2a, the simplified geometry of the neuron-electrode 
interface is depicted. The neuron is modeled as a pillbox 
shaped cell of radius r c =7 um, positioned with an excentricity 
x c , on top of an electrode of radius r e =5 |im. The electrode 
surface is represented by the bottom of a 800 nm deep hole. 
Between the lower membrane surface of the cell and the 
substrate, the sealing gap is modeled (but not over the 
electrode area). 

The geometry of the neuron-electrode interface is 
implemented in ANSYS Finite Element Software release 5.3 
(SAS IP inc., Houston). The medium, surrounding the 
neuron-electrode interface is modeled as a pillbox shaped 
volume conductor of conductivity ov=1.65 S/m, radius 
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fig. 2 (a) Simplified geometry of the neuron-electrode 
interface. The neuron is modeled as a pillbox shaped cell and 
is positioned on top of an electrode with a small sealing gap 
between the cell membrane and the substrate, (b) Finite 
element model of the neuron-electrode interface .The neuron is 
modeled as a cavity inside a volume conductor respresenting 
the surrounding medium. The conductivity of the sealing gap 
is corrected for its actual thickness, (c) The intracellular space 
is assumed to be equipotential and therefore modeled with a 
single node.The membrane is represented by capacitive circuit 
elements, connecting the intracellular node to the nodes on the 
cell upper and lower surface. Finally, the electrode 
impedance and shunt capacity are added to the model using 
circuit elements. 



f m = r t +r c +\x c , positioned at \x c . The neuron is 
modeled as a pillbox shaped cavity inside this volume 
conductor (fig. 2b). The sealing gap could not be meshed 
when implemented with its actual thickness of 10 nm [5]. 
Therefore the conductivity of the gap (the dark ring in fig. 2b) 
was corrected for the difference between modeled and actual 
thickness. 

All volumes in the model are meshed and filled with 
tetrahedral shaped volume elements of type SOLID98 . 
The intracellular potential is assumed constant over the entire 
intracellular space. Therefore, the intracellular potential can 
be represented by a single node. The electrical properties of 
the (passive) neuronal membrane are implemented by 
capacitive circuit elements connecting the neuronal membrane 
to this intracellular node (fig. 2c). All capacities are scaled to 
a value of luF/cm using the area represented by the 
membrane nodes. 

The electrode surface is modeled as an equipotential area by 
coupling all nodal potentials to one principal node in the area. 
A current applied to this principal node will result in an 
equivalent uniform current density over the entire electrode 
surface. The electrode-electrolyte interface impedance is 
added to the model by connection of circuit elements between 
the principal node of the electrode area and a contact node, 
representing U Mt (see fig. lc and fig. 2c). Finally, the shunt 
capacity (between the electrode lead and the medium, across 
the insulation layer) is modeled by a capacitive circuit 
clement. 

The nodes of the upper surface of the volume conductor 
(medium) and of the circular outer surface are set to zero 
potential. A unity current is applied to the contact node. The 
finite element problem is solved at frequencies f= 300 Hz, 
500 Hz, 700 He, 1 kHz, 3 kHz, 5 kHz, 7 kHz, 10 kHz and 20 
kHz, for cell excentricities x c =0 up to x c =4.0 um, in steps of 
0.2 um. The excentricities x c =1.8 um and x c =2.0 urn are left 
out since a proper meshing could not be produced for these 
geometries. 

From every solution the computed potentials at the contact 
node and the intracellular node are saved. Since unity current 
is applied to the contact node, the potential of the contact 
node, U ejlt , equals the total impedance of the neuron-electrode 
system and the potential of the intracellular node, U intra , equals 
the stimulus current transfer function. 



III. Results 

In fig.3 the transition from complete (x c <r c -r e ) to defect 
sealing (x^-r,.) is illustrated by the locations of the 
equipotential surfaces in the volume conductor (sealing gap 
and medium). For convenience, the circuit elements modeling 
the membrane, the electrode impedance and the shunt capacity 
are left out. When the cell is positioned right on top of the 
electrode, all equipotential surfaces are located in the 

scaling gap, which indicates that the potential drops mainly 
over the sealing gap (fig.3a). At x c =1.6 um, the sealing is still 
complete and the potential drop is still located in the sealing 
gap (fig. 3b). When the sealing becomes defect, at x c =2.2 
urn, the potential drop is extended to the rest of the volume 
conductor, indicating a much larger leakage current (fig. 3c). 



The impedance loci (the imaginary part vs. the real part for 
each computed frequency) of the total neuron-electrode 
system are depicted in fig. 4a. The loci are grouped in two 
clusters, which indicate a sharp transition from complete to 
defect sealing. In fig. 4b, some results from previous 
measurements on cultured electrode arrays are redrawn from 
[4] in order to allow comparison with fig. 4a. The measured 
loci are grouped in clusters, similar to the clusters in the 
computed results. This does not surprise, since only a few 
electrodes in the culture were completely covered by a cell. 
The clustering is also present in the computed transfer of the 
applied stimulus current to the intracellular potential, H c (f), as 
can be concluded from the Bode diagram in fig. 4c. For 
stimuli in the frequency range from 100 Hz to 3 kHz the 
transfer function decreases with a factor ten when the sealing 
becomes defect. 

In fig. 5a the sealing resistance is depicted as a function of 
cell excentricity, x c . When the sealing is complete (Xc<2 u.m) 
the sealing resistance decreases only slowly. At the transition 
from complete sealing to defect sealing however, a steep 
descent occurs to a value of less than 100 k£l 
The relation between the stimulus transfer and the cell 
excentricity is depicted in fig 5b for frequencies of 1 kHz and 
10 kHz (these points are also marked in the Bode diagrams of 
fig. 4c). At f=l kHz , this relation is very similar to the 
relation between sealing resistance and cell excentricity. At f= 
10 kHz, the stimulus transfer for complete sealing is lower 
than at f= 1 kHz, but remains more constant until the sealing 
becomes defect. 

From the information presented in fig. 5a and fig. 5b, the 
stimulus transfer can be related to the sealing resistance. 
These results are depicted in fig. 6. The relation appears to be 
linear at f=l kHz. However, at f=10 kHz saturation occurs 
when the sealing resistance becomes higher than 2 MO. 
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fig. 3 Illustration of the transition from complete to defect 
sealing. The colors mark the equipotential surfaces in the 
volume conductor, (a) x c =0. Complete sealing (b) x c =1.8 pm. 
Still complete sealing (c) x c =2.2 fjm. Defect sealing. 




us to intracellular potential 



I 2 

3- 




(a) 

fig. 5 (a) Relation between sealing resistance and cell excentricity. During complete sealing, the sealing resistance remains relatively 
high and decreases only slowly with the excentricity. At the transition to defect sealing, the sealing resistance decreases to less than 100 
k£i. (b)Relation between the stimulus transfer and the cell excentricity at two frequencies. Atf=l kHz, the relation is very similar to the 
relation between sealing resistance and excentricity. At f=10 kHz the stimulus transfer is lower, but remains relatively constant until 
the sealing becomes defect. 
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Although large simplifications are made in the representation 
of the geometry of the neuron-electrode interface, the 
impedance spectra, computed with the finite element model, 
show the same clustering as the measured impedance spectra. 
This clustering is also present in the computed stimulus 
transfer, which indicates that the finite element model 
provides a geometry induced relation between the sealing 
resistance and the stimulus transfer: even a small sealing 
defect annihilates both the sealing resistance and the stimulus 
transfer. 

The highest stimulus transfer that can be reached is likely to 
depend also on the geometry of the neuron, to which large 
simplifications are made. This geometry will modulate the 
relation between sealing resistance and stimulus transfer. 
However, for an accurate geometry based prediction of the 
stimulus transfer, the finite element model can be refined 

x 10 s Absolute transfer of current stimulus to intracellular potential 
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fig. 6 Relation between sealing resistance and stimulus transfer. 
At f=l kHz this relation looks linear. At f=10 kHz the transfer 
esfor sealing resistances higher than 2Mi2 
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V. Conclusions 
The finite element approach offers a powerful method for 
linking the geometry of the neuron-electrode interface to its 
electrical properties. Small sealing defects annihilate both the 
sealing resistance and the stimulus transfer. 
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